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Abstract 

We present a one-dimensional model for diffusion in a fluctuating lattice; that is a lattice which 
can be in two or more states. Transitions between the lattice states are induced by a combination 
of two processes: one periodic deterministic and the other stochastic. We study the dynamics of 
a system of particles moving in that medium, and characterize the problem from different points 
of view: mean first passage time (MFPT), probability of return to a given site (P So ), and the 
total length displacement or number of visited lattice sites (A). We observe a double resonant 
activation-like phenomenon when we plot the MFPT and P so as functions of the intensity of the 
transition rate stochastic component. 

PACS numbers: 46.65. +g, 05.40.FW., 05.10.Ln., 02.50.Eg. 
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A. Introduction 



The problem of diffusion in media submitted to global and/or local fluctuations has 
strongly attracted the attention of researchers. The main motivation was its interest in 
the study and analysis of a large variety of processes. To name just a few examples, we 
consider: random- walks in disordered systems DDQ (ionic conduction in polymeric solid 
electrolytes j^]); transport of Brownian particles that can be in two or more states executing 
a diffusion process in each of them but with different diffusion constants |5l &; resonant 
activation over fluctuating barriers h]j and escape from fluctuating systems [a l2| ; diffusion 
of ligand with stochastic gating 

QLl Q; dynamical trapping problems Qfl HQ- 
Examples of potential applications are charge transport in molecular scale systems 12| and 
photochemical reactions 18{] among others. 

The above indicated phenomena share the property that the switching between the dif- 
ferent configurations or states of the medium is independent of the transport or diffusion 
processes of the walker. Usually, it is assumed that the states are independent of each other 
and that the particles are subject to a Markovian process inducing their motions within each 
state, described by Master or Fokker-Planck equations. 

In previous studies we have considered on one hand the problem of diffusion in fluctuating 
media subject to Markovian and/or non Markovian switchings [3]. On the other hand 
we have studied a stochastic resonance-like phenomenon in a gated trapping system 
Motivated by those previous works, we study here the evolution of particles diffusing in a 
fluctuating medium when the switching mechanism that governs the transition between the 
states has two components: one deterministic and the other stochastic. The parameters 
were chosen in such a way that the effect of the deterministic mechanism alone is not 
enough to produce the switching. However, the simultaneous action of both signals can 
induce the jump over the potential barrier and produce the switching in the medium. Here, 
we consider a periodic signal for the deterministic process, while the stochastic component 
is characterized by a Gaussian white noise of zero mean and intensity £o- I n this way we 
expect to capture some basic aspects of the dynamics of more realistic diffusing systems 
and find situations where the tuning between the stochastic and the deterministic rates can 

induce some resonant phenomena as in Q. 

In order to study this problem we have analyzed the Mean First Passage Time (MFPT) 
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in an infinite system. We have also analyzed the probability that the particle returns to the 
initial position (P so ), characterizing in this way the system's capacity to restore its initial 
condition. We also study a measure of the total length displacement or number of distinct 
visited lattice sites. As it is well known, the MFPT is an extremely important quantity 
because it gives insight into the reaction processes (trapping, etc) and their dependence on 
the details of the underlying dynamics of the systems. The MFPT also provides a useful 
measure of the efficiency of the trapping and has been investigated in a broad range of 
problems from chemistry to biology 20 , 21 , 2, |2| . 



The main goal of this work is the observation of a double resonant activation-like phe- 
nomenon when we analyze the MFPT and the P So as functions of the noise intensity of the 
stochastic component of the transition rate among lattice states. 

The organization of the paper is as follows. In the next Section we describe the model 
in a formal way. After that, we present the results of numerical simulations. In the final 
Section we discuss the results and draw some conclusions. 



B. The Model 

In order to fix ideas we start considering the problem of a particle performing a random 
walk on a continuous or discrete fluctuating medium characterized by N states which are 
labelled by the index j (with j = 1,2,... TV). Following van Kampen 2^| we define the 
probability Uj(t) that the medium has stayed in the j state after a time t since its arrival at 
t = as 

Uj (t) = exp (~J o J2 ^)dt'\ , (1) 

where 7y(£) is the probability per unit time for the medium to jump from level j to level % 
and t is the time it has sojourned in j. 

The switching statistics of the medium %(£); defined as the probability that the medium 
ends its sojourn in the state j after a time between t and t + dt since it arrived at state j at 
t = by jumping to a given state i; is 

v ij (t)dt = u j (t) li] (t)dt. (2) 

It is worth remarking here that Eq. (j2J) is completely general and no extra assumption has 
been made in writing it [24]. As is well known, if the %•(£) are exponential functions of time, 



that is the jij(t) are ^-independent functions (see for instance [20]), we have a Markovian 
switching process between the states of the medium. 

In order to simplify the problem, we assume that the system can fluctuate between 
only two states. The new aspect here is that we have used a combination of noisy and 
deterministic switching mechanisms to describe the transitions between those states. We 
assume that Vij(t) is given by 

v i:j {t) = 9[Bsm{ujt) (3) 

where 6(x) is the step function, and determines when the lattice is in one state or in the other. 
The dynamics of this mechanism is the following: if the signal, composed of a harmonic part 
plus the noise contribution £, reaches a threshold £ c the lattice changes its state, otherwise it 
doesn't. We are interested in the case where £ c > B, that is, without noise the deterministic 
signal is not able to induce a change of the lattice s state. 

Finally, in order to complete the model, we must give the statistical properties of the 
noise £. We assume that £ is an uncorrelated Gaussian noise of width £ > i- e - 

(£(*)£(«)) = Mo- (4) 



Note that this is not a standard white noise j2l| due to the Kronecker symbol (instead 
of Dirac ' s delta 5(t — u)). For each fixed state j of the medium the transport process is 



Markovian and we denote its corresponding "propagator" [24 1 by Aj. These propagators 
are differential operators (matrices) in the case of a continuous (discrete) medium and its 
structure depends strongly on the character of the fluctuations. 

In we have been able to obtain some analytical results for the MFPT in the case of 
Markovian and non-Markovian processes that were satisfactorily compared with numerical 
simulations. However, and due to the complex dynamics of the present case, we were not 
able to solve this problem analytically and have resorted to MonteCarlo simulations. 

The numerical simulations were performed on a one dimensional infinite lattice. The 
particles were initially located at the lattice site sq in a given state. We assumed that the 
propagators Aj were arrays W describing a discrete one dimensional random walk in each 
state of the lattice, and also assumed jumps only to first neighbors. The jumps within each 
state were characterized by two parameters: Aj, the temporal rate of jump and rjj, the "bias" 
to make a jump in a given direction. For the case in which we studied the return to the 



origin, the system started its motion at the origin. All simulations shown in the following 
figures correspond to averages over 100000 realizations. 

In the next Section we present several figures where we show the results of our simulations. 
In the whole study we have used state 1 as the standard state, characterized by the following 
parameters: Ai = 1 and 771 = 1 (we have chosen the bias to the left side of the lattice), while 
state 2 is the test state characterized by A2 = l.We have chosen the same temporal rates 
because we want to focus on the effect of the bias and the noise intensity on the system 
response. Hence, rj 2 together with £ , are the parameters we varied. In addition, and for 
all figures, we have adopted the following parameters for the deterministic signal: frequency 
u) — 1., and amplitude 5=1.; while the activation threshold is £ c = 2. 

C. MonteCarlo Results 

In Fig. Q we show the results for the MFPT to reach the origin (s = 0) for a particle 
that started its motion at site so = 4, as a function of the noise intensity £o- From this 
figure, it is apparent that there is an optimal value £0 f° r which the MFPT has a minimum 
(or viceversa, the "activation rate" has a maximum). This behavior, that is a manifestation 



work. When £0 is small, the transition rate is low, and the particle is in an "unfavorable" 
state, it will remain there for a long time, contributing to a larger MFPT. When £0 is large, 
the transition rate will be high, and the particle will not show a net motion, again yielding 
a large MFPT. But there is an optimal £0 where the combination between transition rate 
and motion yields a minimum of the MFPT. It is worth remarking here that this situation 
occurs when the bias in each state is strong enough but they are in opposite directions, 
however the phenomenon disappears when the tendencies in both states are similar. 

Figure El shows the dependence of the MFPT on 772 for two values of £0 corresponding 
to low and high noise intensities. In this figure we can recognize two main aspects. Firstly, 
the MFPT decreases as r/ 2 grows, which is a logical result because increasing 772 means that 
we have an increasing bias to the origin. Secondly the figure shows an important difference 
in the MFPT for small r]% , where we see that the noise intensity plays an important role. 
For small £0 the transition rate between the states is low and the particle can remain a long 
time moving to or away from the origin. For large 772, the differences between the behavior 



of the so called resonant activation phenomenon 
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FIG. 1: MFPT vs £o> f° r % = 0.1(open circles) and rj2 = 0.05 (black-squares). 

for different noise intensities disappear. Also, at intermediate values, the MFPT remains 
higher for a smaller noise intensity. 

Figure El depicts P SQ , the probability of return to the origin, as a function of rj2- The 
initial condition of the system was at site zero (sq = 0) in state one, and the simulation time 
was 10000. For comparison, the case of only one lattice is also shown in this figure. The 
latter presents a maximum for an intermediate value of r\2 > that is when there is no privileged 
direction, otherwise the probability decreases as the particle is forced to move away from the 
origin. For the two state case, we can see that, for small noise intensities, the noise enhances 
the system response for a fixed (small) value of r\2- This is due to the possibility of changing 
from the unfavorable state to the one where the particle has a larger probability of moving 



6 



1000 ■ 



"i 1 1 1 1 1 1 1 r 



O, 



O 



o 



100 - 



o 



o 



_Q 

03 



o 



o 



o 



o 



10 - 



0. 



® 



— I ' 1 ' 1 ' 1 1 1 ' — 

0.0 0.2 0.4 0.6 0.8 1.0 



FIG. 2: MFPT vs 772 , the bias in state 2. The circles corresponds to the MFPT values for £0 = 1-0, 
and the squares for £0 = 10.0. 

towards the origin. If we increase the noise intensity, the response is enhanced, that is the 
return probability becomes larger. However, if we continue increasing the noise intensity 
finally the response for low r/2 is reduced and there appears a maximum in the probability 
indicating that there is an optimal value of 772- For large 772 values, the response becomes 
smaller for all noise intensities as expected. Here we have another important result of our 
study, a kind of double resonant behavior both in £0 and 772. 

Figure |U depicts P so as a function of the noise intensity £ , for different values of r] 2 . The 
increase of this probability with increasing £ is apparent, reaching a kind of plateau for 
£0 > 10. However, for rj 2 = 0.3, a maximum around £0 = 10 is insinuated, resembling a 
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FIG. 3: P SQ vs 772 - The circles represents the one-state case, the squares are the data for £0 = 5.0, 
the crosses for £0 = 10.0 and the triangles for £0 = 20.0. 

resonant like phenomenon. The fact that this probability decreases with increasing 77 2 is 
again a logical result because larger r/ 2 implies an increasing bias to the origin. 

Figure 03 shows the total length displacement or number of different lattice sites visited, 
that we indicate with A, as a function of the noise intensity £0, f° r different values of 772. For 
small noise intensity there are no differences among different values of 7/ 2 , but for increasing 
£0 the number of visited sites is strongly reduced for 77 2 = 0.3, 0.5, but remains still high for 
r] 2 = 0.9. This is again consistent with a larger bias to the origin with rj 2 larger. 

Figures El and |H1 show the total length displacement or number of visited lattice sites, 
indicated with A, as a function of r/ 2 and for different values of £o- In all three cases the 
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FIG. 4: P So vs £o- The circles represents the data for 772 = 0.3, the squares for 772 = 0.5 and the 
crosses for 772 = 0.9 

simulation time was t = 1000, and the initial state of the system was the state 1. In Fig. |U] 
we compare different cases where the noise intensities differ by an order of magnitude. When 
£0 — 1 (more precisely £0 = 1.5), the number of visited sites remains almost constant with 
772. For £0 = 10, the curve starts at a small number of sites and increases monotonically, 
almost linearly, with 772- For larger values of £0 (£0 = 100), the presence of a minimum is 
apparent. A remarkable fact is that before and after the minimum in A, A behaves linearly 
with r/ 2 . It seems that the high number of visited sites for small £0 when r/ 2 is also small, 
reduces with increasing £ before rising again for still larger values of £o- This argument was 
supported by complementary simulations, see for instance Fig. [7| where we have depicted 
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FIG. 5: Ao vs £o- The simulation time in this and the following figures was t = 1000. The circles 
represents the data for 7/2 = 0.3, the squares for 7/2 = 0.5 and the crosses for 772 = 0.9 

the same situation but now for a lattice with only one state, as well as for the two-state 
lattice and in the range of large noise intensities (£ = 20,30, 100). The indicated effect is 
apparent. For the one-state lattice, the minimum is located at i] = 0.5. The reason is that 
i] = 0.5 corresponds to the situation of unbiased particle motion, that is we have the same 
probability of moving to the left and to the right of the lattice. This is a point of symmetry 
in the sense that a larger or smaller 77 implies a biased motion in one direction or the other 
(changing rj by 1 — i] only changes the direction of motion but not the number of visited 
sites). 

We can see how the noise changes the position of the minimum in the two-state lattice 



10 



1000 - 



900 



800 - 



700 - 



600 



< 500 - 



400 - 



300 - 



200 - 



100 



~\ ' 1 1 1 ' 1 1 1 ' r 

++++++++++++++++++ 



■O 



■o 



o 



o 



o 



o 



o 



o 



o 



o 



o 



o 



o 







o 



o 



cP 



-i 1 1 1 1 1 1 1 1 1 1 1 r- 

-0.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 



FIG. 6: A vs r/2. The crosses represents the data for £o = 1.5, the squares for £o = 10.0 and the 
circles for £o = 100.0. 

relative to the one-state case. The strong dependence of the system ' s response for small 
bias is apparent, as well as the linear scaling of A with 772 before and after the minimum . 
The decrease in the value for small rj 2 is also shown in Fig. where we have depicted the 
same situation as before but now for £0 = 1.5,3. and 5.. The indicated decrease is clearly 
seen. From this figure we can see that the indicated minimum disappears for small noise 
intensities, and the "full" linear dependence of A on 772 in this range. 
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FIG. 7: A vs 772- The result for a one-state lattice is indicated by circles. The black-squares 
represents the data for £0 = 100.0, the black-triangles for £0 = 30.0 and the crosses for £0 = 20.0. 

D. Conclusions 

We have studied here the evolution of particles diffusing in a fluctuating medium when the 
switching mechanism that governs the transition between the states has two components: one 
deterministic and the other stochastic. We have chosen a periodic signal for the deterministic 
part while the stochastic component is characterized by a Gaussian white noise of zero mean. 
The parameters were chosen in such a way that the effect of the deterministic mechanism 
alone is not enough to produce the switching, but the simultaneous action of the deterministic 
and the stochastic signals can induce the jump over the potential barrier and produce the 
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FIG. 8: A vs r/2- The circles represents the data for £o = 1-5, the squares for £o = 3.0 and the 
crosses for £o = 5.0. 

switching. The problem was characterized studying the MFPT in an infinite system, and 
we have also analyzed the probability that the particle returns to the initial position, and 
the total length displacement or number of visited lattice sites. 

The main goal of this work is the observation (see fig. ^) of a resonant activation-like 
phenomenon when considering the MFPT as a function of the noise intensity of the stochastic 
component of the transition rate. Another important result is the kind of double resonant 
behavior, both in £ and t? 2 , observed in fig. El 

It is also worth remarking here the behavior of the number of visited lattice sites (A) as 
a function of r] 2 , with £ as a parameter. For the case of the one-state lattice, a minimum 
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in the number of visited lattice sites exists for rj = 0.5. For the two-state lattice, that 
minimum changes with the noise intensity. For large values of £o> it shifts but is still well 
marked. However, when £o is reduced the minimum is less marked and finally, for very small 
values of £o, it disappears. This is another example of the "constructive" role that noise can 
play, in this case associated to diffusion on a lattice with several (here two) states. Another 
remarkable fact is the "piece-wise" linear dependence of A on 772 

According to the discussion above on the motivations of the present work, a natural step 



in this research of diffusion in a fluctuating medium will be, as in Ref. |l6j |. the inclusion 
of a gated trap mechanism at the origin. The aim will be to analyze the interplay between 
both deterministic-plus-stochastic mechanisms, the switching and the gating one, and the 
coherent or interference effects that could arise. This problem will be the subject of further 
work. 
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